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Abstract 

We show that the determinant of a Hankel matrix of odd dimension n whose entries are the 
enumerators of the Jacobi symbols which depend on the row and the column indices vanishes iff 
n is composite. If the dimension is a prime p, then the determinant evaluates to a polynomial of 
degree p—1 which is the product of a power of p and the generating polynomial of the partial sums 
of Legendre symbols. The sign of the determinant is determined by the quadratic character of —1 
modulo p. 

The proof of the evaluation makes use of elementary properties of Legendre symbols, quadratic 
Gauss sums and orthogonality of trigonometric functions. 
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1 Introduction 

For an odd integer n, and fc = 1, 2, . . . , n define the polynomials 

k 



flfc 



(x) — J(k — m, n)a 



in which J(a, m) is the Jacobi symbol defined for odd integers m by 

\piJ \p2^ \pk-' 

where the prime factorization of m is m = ' ' 'Pit ^^"^ ^ prime p, is the Legendre symbol 

defined by 

r ifp|a, 
— ^ = < 1 if a is a quadratic residue mod p, 

I —1 if a is a quadratic nonresidue mod p . 
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For example when n — 3, the first five polynomials are 



ai{x 
a2{x 
as,{x 
04(2; 



= 1 
= x-1 



x^ -x'^ + 1 

x'^ - x^ + x - 1 



It is easy to see that ak{x) is a monic polynomial of degree k — 1 and 0^,(0) ~ J{k,n). Consider the 
n X n Hankel determinant 

Hn{x) = det[ai+j_i(a;)]i<ij<„ . 



(1) 



As an example, 



Hs^x) = dct 



1 



x-1 



X — 1 X 



X — X 
„3 ^2 



X^ ~ X — .T^ + 1 X^ — 2 

A few other determinant evaluations for small n are as follows 



x'^ — X +1 

+ a; - 1 



H^x 
Hj{x 
H^ix 

Hi3{x 
Hi5{x 
Hn{x 

HMx) 



5x'^{x-l){x + l) 

-A9x^{x^ + 2x^ + x^ + 2x + l) 



-14641a;2(a;^ + x^ + 2x^ + 3x^ + 2x^ + x^ + 1) 

371293a;^(x - l)(x + l)(a;^ + 2x^ + 2x^ + 3a;'* + 2a;^ + 23;^ + 1) 



410338673a;2(a;-l)(x + l) 

{x^^ + 2a;" + 2x^° + Ax^ + 3x^ + ix^ + 2x^ + Ax^ + 3x'^ + Ax^ + 2x^ + 2x + 1) 
-16983563041a;2 



(x 



16 



+ a;^2 + 2a;" + 3x*° + 2x^ + 3x* + 2x^ + 3x^ + 2x'^ + a;"* - a;^ + 1) 



Recently, Chapman [2] evaluated Hankel determinants of certain 



dimensional 0-1 matrices 



built up from the Legendre symbol defined modulo a prime p. These evaluations give 



det 



= det 



i<«j<- 



1 - 



= - 1 



for any prime p > 3, p = 3 (mod 4). [2] also includes additional conjectures related to such determi- 
nants. In this paper, we prove the following evaluation of the Hankel determinant Hn{x): 

Theorem 1 Hn{x) identically vanishes unless n ~ p is a prime. For p prime, 

p-i 

Hp{x) = (-l)^p'^^6fca;*'- 

fc=0 
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where 




) 



(2) 



Furthermore, Hp(x) is divisible in ^[x] by forp = 3 (mod 4) and by x'^(x^ — 1) for p — 1 (mod 4). 

The properties of the Jacobi and Legendre symbols and Gauss sums that we make use of in the 
proof of Theorem[T]can readily be found in most books on number theory: we mention only [T], [5], [3]. 



We divide the proof of the theorem into a series of lemmas, and start with recording the following trivial 
property of the polynomials ak{x): 

Lemma 1 



2.1 The composite case 

Now we show that Hn{x) = iff n is composite, and then determine the structure of Hp{x) for p prime. 
Lemma 2 Hn {x) identically vanishes for n composite. 

Proof Let r,; = (0^,0^+1, . . . ,ai+„_i) denote the i-th row of the matrix in ([1]). Let e^ denote the 
n-dimensional unit row vector with 1 in the i-th coordinate and elsewhere, with e* denoting its trans- 
pose. The proof is in two cases depending on whether or not n is a perfect square: 

Case I: n = is a perfect square. 

We claim that in this case the four rows ri, r2, fm+i, rm+2 are linearly dependent. More precisely 



2 The proof of Theorem [T] 



<ik+i{x) = Jik + 1, n) + xa^ix) . 



V2 - XVi = r,„+2 - XTm+l ■ 



From Lemma [U 



2 



m' 




(3) 



and 



2 




(4) 



Note that 




gcd{a, to) > 1, 
gcd{a, to) = 1 . 



Since 



gcd{i + 1, m) = gcd{i + m + 1, to) , 
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the right hand sides of ^ and (jH) evaluate to the identical 0-1 vector. 



Case II: n = p'^'^'^^q with p prime, p Jl q. 



Let m — p^'^^^ . In this case we show that the following linear dependence among the rows holds: 

p-i 

^(r,,+2 - XTiq+i) = . 
1=0 

By Lemma [T] the j-ih. entry of the vector on the left is 



j=0 



^J{iq + i + l,mq) = .J(j + 1, q)^ J{iq + j + l,m) 

1=0 

p-i 

= J(j + l,g)^ J(2(z + j + l,p) 



2 = 

^0 + 1''?)E(^)-0 



2.2 The prime case 

Let now n = p be prime. Using Lemma[T]and replacing r^+i by r^+i — xvi for i = 1, 2, . . . ,p, we obtain 

ai{x) a2{x) ■■■ ap^i{x) ap{x) 



Hp{x) — det 



(f) (I) 
(f) (i) 



p 



(^) (^) 



(f) (^) (^) (^) 

Since afe(x) is of degree fc — 1, Hp{x) is a polynomial of degree p — 1. 

Consider the p x p matrix 

(1) © 
- (^) 



A — 



- - 



(l) (Ir^) ■■• (^) 



j + i - 1 



(5) 



(6) 



Let denote the cofactor of the entry («, j) of Ap. Expanding the determinant in ([5]) by the first 
row, we have 

p 

Hp{x) = ^cijaj(a:) 

3 = 1 
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and the coefficient of the leading term is the cofactor ci.p = det Cp where Cp is the (p — f ) x (p — f ) 
matrix 

(?) 



(3) (f) 



p 



P 



(7) 



l<'i,j <p— 1 



(f) (^) ••• 

First we show that the ci.j's, and in fact all cofactors of Ap are identical. 
Lemma 3 All cofactors of the matrix Ap are identical. 
Proof We note that Ap is a symmetric matrix with the i-th row sum 

'i+j-V 



E 



p 



= 



for every i. Since the row sums vanish, the cofactor c^j is independent of j. By symmetry, Cij is also 
independent of i. One way to prove Lemma [3] combinatorially is to use the standard weighted version 
of Kirchoff's matrix-tree theorem [7], [8j. We include it here for completeness. Consider the complete 
graph Kp on vertices {1, 2, . . . and introduce the indcterminates Xij for 1 < i, j < p- Define 



and define the weighted Laplacian matrix by setting Lp = [Lij]i<ij<p, with 



— 



D, 



if i 

if i^ j 



Let Sp {Kp) denote the set of spanning trees of Kp. For any given index i, we can consider a.T e Sp {Kp) 
as being rooted at vertex i. This simply gives an orientation to each edge e = {r, s} of T by orienting 
it from r to s iff s is closer to the root than r in T. Define the weight of e € T by Wi{e) = Xr,s and the 
weight of T itself by 



eeT 



Then any cofactor Cij of an element in the i-th row of Lp is identical and evaluates to 



TeSp{Kp) 



(8) 



This is the content of the weighted generalization of Kirchoff's matrix-tree theorem. Suppose we 
specialize each Xr,s, r ^ s to a numerical value such that the resulting matrix is symmetric (i.e Xr,s 
and Xs,r are assigned the same value). Given a T G Sp{Kp), Wi{T) then specializes to a fixed value 
independent of i since the symmetry of the matrix implies that either edge orientation results in the same 
numerical weight for the edge. Therefore the sum in ([5]) evaluates to the same quantity independently 
oii,j. • 
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Since ci.j = det Cp for all j, we have proved 
Lemma 4 

p 

Hp{x) — det Cp Qj (x) 

Note that 

p p-i 

i=l /c=0 

where is as given in Next we evaluate det Cp. 



Lemma 5 



detCp = (-l)V^pV^ . (9) 



Proof Let Ep denote the p x p exchange matrix which has 1 's along the anti-diagonal and O's elsewhere. 
Clearly, 



detEp = (-1)^^^^ 



Let Bp = CpEp. Then 



Bp 



l<i,j<p— 1 



and 

detCp = (-l)T^ detBp . (10) 

Note that Bp is symmetric for p = 1 (mod 4) and skew-symmetric for p = 3 (mod 4). 

We determine the spectrum of Bp, and compute det Bp as the product of its eigenvalues. This 
results in the evaluation of det Cp that we need through (jlOp . 

Let / — V — 1 and ^ = e p denote a primitive p-th root of unity. For 1 < r < p — 1 consider the 
Gauss sum 

Then 

,.. = ( (?)^ (U) 

[ r{i)VP a ps3 (mod 4) 

A proof of Gauss's evaluation of gr can be found in 6J. Changing the summation index, we can write 

9r-f.C^)c^'-''- (12) 

We will give the details of the proof for primes of the form p = I (mod 4) . The proof for primes p = 3 
(mod 4) is similar. 

For p= I (mod 4), is even, and det Cp = det Bp. Using pT|) and (fT^ . we have 

£ f ^) <-" - (;) 
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or 



E ( V ) c-^^ = ( ) vpc- . (13) 



Equating the imaginary parts in (|13p . 

sm = - sm 

\ P J P \PJ P 

Therefore for every r which is not zero modufo p, the vector 

p-i 



. 2nrj \ f 
sm e' 



is an eigenvector of Bp corresponding to eigenvalue j The vectors corresponding to r and p — r 
differ only in sign. Therefore if we let 

Ti = {ur\l<r< 

then exactly half of the Ur € Ti are eigenvectors of Bp corresponding to eigenvalue y/p, and the other half 
are eigenvectors corresponding to the eigenvalue —^/p- For 1 < r < s < we have the trigonometric 
identity 

p— 1 f 

E27rrj 2Trsj if r < s 

sm sm — ■! 

P P [ 2 '^^"^ = ^ 

where the r = s evaluation is a consequence of the general trigonometric identity 

E. 9 n cosm + 1 kc sm 
snr = ^ „ / (14) 
2 2sinx ^ ' 

O) P- 30). Therefore the 2zl eigenvectors in Ti are orthogonal, and so linearly independent. 

Next we obtain a set of — 2 more eigenvectors of Bp. Equating the real parts in p^. we obtain 

Let 



^ ^ ^ 27rrA f 
cos e ; 



These are not themselves eigenvectors because of the extra term in (|15p . But the nonzero vectors 
of the form 

Vr - Vs (16) 

for l<r<s<p— 1 are eigenvectors of Bp as long as = {j;^ ■ We will single out 

^-2 
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of these eigenvectors, half corresponding to the eigenvalue y^, and the other half to — y^. Let 5 be a 
generator of the multiplicative group 2*. Then = —1. Put h = g^. Since 

taking r = 1, the — 1 vectors 

for k — 1,2, ... , — 1 are eigenvectors of Bp corresponding to the eigenvalue y/p. Similarly, 
and taking r = g in (fTB|) . the — 1 vectors 



PJ \ P 



Vg - Vghk 

for fc = 1, 2, . . . , 2^ ~ 1 are eigenvectors of Bp corresponding to the eigenvalue — y^. These eigenvectors 
are of the form 

Vl — = \ cos COS 1 



in the first case, and 



in the second. Let 



p-i 



27rgj 27r5ft,'''j 



— Vghk = I COS — — — COS — - — - I e^- 



p — 1 p — 1 
T2 = {vi - Vhk I fc = 1, 2, . . . , — 1} U {vg - Vg^k I fc = 1, 2, . . . , — 1} . 

Finally, consider the two vectors 

Thus wi is a 0-1 vector with a 1 for every index for which the row sum of Bp is 1. Similarly, W2 is a 
0-1 vector with a 1 for every index for which the row sum of Bp is —1. 

The fact that wi is an eigenvalue of Cp (and also of Bp) is a consequence of the identity 

p-i 



E 
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To prove this identity, write it in the form 
In this latter form the identity can be proved by expanding the left hand side and making use of 

3-1 

' t+J \ (J 



which holds for p ^ i from the general orthogonality condition 

p— 1/. 7\/- ;\ ( 1 -r- 

i + k \ f J + k\ I p — I it I = J 



For Bp, we obtain 



fe=0 



Ei^K^iH":; ;;;;; m 



so that wi is an eigenvector of Bp corresponding to eigenvalue 1. Similarly, W2 is an eigenvector of Bp 
corresponding to eigenvalue —1. Putting 

T3 = {W1,W2} 

we have p — 1 eigenvectors in 

T1UT2U n 

with corresponding to eigenvalue y/p, corresponding to eigenvalue ~y/p, and one each for the 
eigenvalues ±1. To show that there is no linear dependence among these vectors, we proceed to show 
that any two vectors u,v E T1UT2UT3 are orthogonal. We have already done this for u,v E Ti. For 
u,v E T2, we need to show 

p-i 



y I cos cos 1 I cos cos 1 = 



p-1 

E 

p-1 



p p J \ p p 



2T:gj 2TTgh^'j\ ( 2T:gj 2'Kgh'^j 
cos cos 



p P J \ P P 



^ , 2nj 2Tdfj\ ( 2TTgj 2TTgh'j \ 
y cos cos cos COS = 

~{\ p p J \ p p J 

for r ^ s (mod p). 

These identities follow from 

n-i I ^1 iis^r,n — r 

E2Trrj 2iisj I 

cos cos — < n— 1 ifs = r = § 

n n \ „ ^ 

J=i { ifs = l,n-r 

which holds for l<r<s<n— 1, and generalizes the twin identity to ([1 



n 

COS ]x = — h ij cos nx sin(n + Ijx esc x 



3 = i 
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m,P- 31). 

To prove that the vectors in Ti are orthogonal to the vectors in T2, we use the orthogonaUty relations 

E" ^ 2TTrj 27:3] 
cos sin = 

n n 

valid for all integral r,s,n. 

Finally, below are the identities that are needed to prove that the vectors in T3 are orthogonal to 
vectors in Ti and T2. If p is a prime of the form 4A; + 1, then 

E «in^=0 (18) 

i = i 
-) = 1 



for any r, and 



E = (19) 



(i) = 1 



and 



The first one of these can be written as 



Clearly, 

sm = 

P 

i=i ^ 

by looking at sin as the imaginary part of C and summing the geometric series in C^. Therefore to prove 
p^ . it is enough to prove 

y(l\ sin^=0 
^ VP/ P 

which is an immediate consequence of the evaluation of the Gauss sum by equating the imaginary parts. 
The identities ^9]) and ((20)) are obtained by evaluating 

again by making use of the evaluation of Gauss sums. We sum the geometric series in and equate the 
real parts. 
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Therefore the spectrum of Bp consists of ±l,±^/p where 1 and —1 each have muhiphcity one, and 
yj^ and —\/v each have muhiphcity — 1. This gives 

detBp = {-Vf^-p^ . 

This completes the proof of Lemma [S] • 

I am grateful to the anonymous referee who suggested an alternate, and somewhat more economical 
proof of Lemma [5l I would like to sketch this approach here. Let 

-3 



Dp 

We can view Bp as a submatrix of of Dp obtained by deleting the first row and column of Dp. Dp is 
a circulant matrix, and therefore it has a basis of eigenvectors consisting of the (1, C", C^"^ , . . . , (^(p^^)''). 
The eigenvalues are (p— l)/2 of them equal gi, (p— l)/2 of them equal —g\ and also go = must be 
included. If we have an eigenvector of Dp with first entry zero, deleting that zero gives an eigenvector of 
Bp with the same eigenvalue. Taking differences of the above basis elements gives (p — 3)/2 independent 
eigenvectors of Bp with eigenvalue gi, and (p — 3)/2 with eigenvalue —g\. This accounts for all but two 
eigenvectors of Bp, and these two are w\ and w^- 

Remark 

For p = 3 (mod 4), the spectrum of Bp consists of ±/, ily/p where / and —I each have multiplicity 
one, and ly/p and —ly/p each have multiplicity — 1. In this case Bp is skew-symmetric, so the 
determinant is non-negative. 



3 Special values 

We can obtain factors of Hp{x) by finding zeros of 



(21) 

k=0 

where bk is as given in 

Lemma 6 For any p, \ Hp{x). 

Proof It is easy to see that for any odd prime, ~ bi = 0. Therefore Hp{x) is divisible by x'^. • 

Next we consider the case p = 1 (mod 4) . 
Lemma 7 If p = I (mod 4), then we also have (x^ — 1) | Hp{x). 

Proof The polynomial (pij) evaluated at a; = 1 and a; = — 1 are 

P ^ / 1 



k=l m=0 
p k 



P 



fc=l m=0 
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respectively. We will show that both of these evaluate to 0. Rearranging the first sum, 

p k / _ \ 

EE (^) = E(p-^"+i) 

k=l rr, 

Therefore it suffices to show that 



P J — \ P 

k=l m=0 \ ^ / m=l ^ ^ 



E - (7) - 0' (22) 

m=l \ ^ / 

i.e. the sum of the quadratic residues minus the sum of the quadratic nonresidues mod p vanishes. 
Since 

the map m ^ p — m permutes the quadratic residues among themselves, and the nonresidues among 
themselves. Since each of these sets have an even number of elements for p = 1 (mod 4), this map has 
no fixed points. Therefore both the residues and the nonresidues mod p sum to 

and (|22p follows. The second sum in can be rearranged as 

p— 1 p-1 

' 2m — 1 \ , , p^-i \ ^ / m 



E F^ =-(-i)"^E - (23) 



p J , \ p 

and in this case the map m ^ p — m shows that there are equally many residues mod p in the range 
{1, 2, . . . , ^y^} as in the range {^^ + 1, . . . ,p— 1}. A similar statement holds for nonresidues. Therefore 
the right hand side of ((23l) is zero and Hp{x) is divisible by a;^ — 1 for p = I (mod 4). • 

Note that the elementary arguments we gave for the proof of the evaluations in Lemma[7]can directly 
be obtained from the following result (see [9], also [3]): 

Proposition 1 Let p be an odd prime and suppose F is a complex-valued function defined on the 
integers, which is periodic with period p. Then 

p-1 , 

J 



E^(^)+EUr(^-) = E^(^")- 

Finally, we remark that the coefficients of the quotient polynomials 

(-l)^Hp(x) (-l)^iJp(x) 



p 2 X'' p 2 x'^yx'^ — 1) 

over '2.[x\ can be written in terms of the partial sums hk- For p = 3 (mod 4) these coefficients are 
simply bk+2- For p = 1 (mod 4) the coefficients are partial sums of odd or even indexed hi, depending 
on the parity of k. 
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